Surface segregation of conformationally asymmetric polymer blends 
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We have generalized the Edwards' method of collective description of dense polymer systems in 
terms of effective potentials to polymer blends in the presence of a surface. With this method 
we have studied conformationally asymmetric athermic polymer blends in the presence of a hard 
wall to the first order in effective potentials. For polymers with the same gyration radius Rg but 
different statistical segment lengths Ia and Ib the excess concentration of stiffer polymers at the 
surface is derived as 8pA{z = 0) ~ {l~^ — l~i^)\n{Rg /l^), where Ic is a local length below of which 
the incompressibility of the polymer blend is violated. For polymer blends difi'ering only in degrees 
of polymerization the shorter polymer enriches the wall. 

PACS numbers: 61.25.Hq, 68.47.Pe, 83.80.Tc. 



I. INTRODUCTION 

The effects of surfaces on the behavior of polymer melts 
and blends is of basic importance in their numerous ap- 
plications such as adhesion, lubrication, wetting, cataly- 
sis, etc. Q. The structure and properties of the blends 
and other polymeric materials within a few nanometers 
at a surface can differ significantly from corresponding 
properties in the bulk. For example, in polymer blend 
a segregation of one of the components to the surface is 
possible even if the blend is miscible in bulk. Hence the 
questions how the surface can change the properties of 
polymeric materials and how it may be controlled are of 
practical interest. Despite the large theoretical and ex- 
perimental interest in the behavior of polymer blends in 
the presence of surfaces and achieved basic understanding 
there is no satisfactory analytical treatment of segrega- 
tion of polymers at surfaces. Most studies of polymer 
blends near surfaces are based on phenomenological ex- 
pressions for a free energy, which include surface terms 
that account for adsorption or repulsion of a particular 
type of monomers 2, 3]. Minimization of the free energy 
gives equilibrium concentration profiles for each compo- 
nent. There exist more rigorous approaches, which allow 
one to derive the concentration profiles starting from the 
microscopic polymer statistics in the presence of a sur- 
face. One of such microscopic approaches is the integral 
equations method which can be applied to various site- 
site or hard-particle models of a dense polymeric system 
0. This method having many advantages such as abil- 
ity to predict microscopic correlations between different 
types of monomers and between monomers and surfaces 
requires a considerable amount of numerical computa- 
tions. Most of analytically treatable methods, which rely 
on the continuum Gaussian chain model, take into ac- 
count the monomer-monomer interactions usually using 
either the random phase approximation |5| j which is most 
suited to treat systems in the weak segregation limit or 
self-consistent-field theories 0, 0| , which are most suited 
to treat systems in the strong segregation limit. 

Recently there has been a big deal of attention on the 
surface segregation due to conformational asymmetry of 



the molecules of the polymer blend and due to differences 
in topology It was 

established that the composition of polymer blends in 
the vicinity of surfaces can be different from the bulk 
composition even for neutral surfaces. It was found that 
for polymer blends composed of polymers with differ- 
ent degree of polymerization but of chemically identical 
monomers, shorter polymers are in excess at the wall. 
It was also demonstrated in simulations 0, lllj and sup- 
ported by calculations using the integral equation the- 
ory |0| that stiffer polymers are present in excess in the 
vicinity of the wall. However the self-consistent field the- 
ory developed in Ref. predicts the opposite effect, i.e., 
the excess of more flexible polymers. Unfortunately, no 
predictions on the behavior of polymer blends of chem- 
ically identical polymers with different degrees of poly- 
merization were made in Ref. [T^ . 

In this paper we present the analytical study of the 
behavior of athermic polymer blend in the presence of a 
hard wall using the generalization of the Edwards' col- 
lective description of dense polymer systems in terms of 
effective potentials Tj, "20] to polymer blends in the pres- 
ence of a neutral surface. The bare one-polymer Green's 
function G obeys the Dirichlet boundary condition. We 
show that a partial summation of graphs results in replac- 
ing the bare G with the effective one-polymer Green's 
function Gr^ which, as we argue, obeys the reflecting 
boundary condition. The bare and effective Green's func- 
tions are related by the Dyson equation, where the self- 
energy E is defined by series of graphs. This part of 
our work is similar to Ref. ^3 ) however with the follow- 
ing significant difference. In the present work the Dyson 
equation results in an integral equation for Gr , which de- 
termines the relevant reference state to describe polymer 
melts in the presence of a neutral wall. The concentration 
profiles are due to fluctuations, which are not taken into 
account in self-consistent field theories. The method we 
use can be applied in a straightforward way to study the 
behavior of polymer blends and copolymer melt in the 
presence of selective surfaces, the dimensions of polymer 
molecules in the melt, the distribution of polymer ends, 
etc. 
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The paper is organized as follows. Section lll Al outlines 
the statistics of a single polymer chain in the presence of 
a hard wall. Section FlI Bl introduces to the collective de- 
scription of dense polymer system. Section lll CI contains 
the discussion of the behavior of the effective potentials 
and screening effects in the presence of a hard wall. Sub- 
section IIIDI introduces to the collective description in 
the presence of a neutral surface. Section IIIII contains 
calculations of the excess monomer concentration of con- 
stituents of an incompressible athermic polymer blend 
in the vicinity of a hard wall. Section IIVI contains our 
conclusions. 



II. COLLECTIVE DESCRIPTION OF DENSE 
POLYMER SYSTEMS 

A. Polymer chains in the presence of a hard wall 



given by 

G(r,iV;r') = Goir,, ~ r^, N;0) [Gq{z - z' , N;Q) 

+ Go{z + z\N;0)]. (6) 

The Laplace transform of the z part of the Green's func- 
tion with respect to N, which we will need in the follow- 
ing, is given by 



Go{z - z',p) = 



1 



2a^ 



exp 



z~z'\^). (7) 
a 



The monomer density of a single polymer chain 

n{z,N) = / ds{S{z - z{s))) 
Jo 



(8) 



can be expressed through the Green's function of the 
polymer chain as follows 



The Green's function of a free polymer, which is pro- 
portional to the relative number of configurations of the 
ideal chain with the ends fixed at r and r', and gives un- 
der appropriate normalization the distribution function 
of the end-to-end distance, obeys the Schrodinger type 
differential equation 19i..21i| 



d 
dN 



2v72 

a Vr 



G{r,N;r') = S{r-r')S{N), (1) 



where N is the number of statistical segments, and 
— P/6 with I being the statistical segment length 
of the chain. The distribution function of the end-to-end 
distance obtained from Eq. reads 



N 



oo oo 



n{z,N) = J ds J dz' J dz"G{z',z,N ^ s)G{z,z",s). 

" 

(9) 

The straightforward computation using the Green's func- 
tion obeying adsorbing and reflecting boundary condition 
yields 

n{z,N) = [2erf(z/2) + z2[l + erf(z/2)] 

2z o 
-I- exp(-z'^/4) - erf(z) 



2z 

2z^erf (z) -= cxp(-z^) 



(10) 



Go(r,iV;0) 



/ 1 



3/2 



exp 



(2) 



In the presence of a hard wall we have to impose an 
appropriate boundary condition on Eq. (QJ . For polymers 
in a dilute solution with a hard wall situated at z = 
one should use the Dirichlet boundary condition 



G(r,iV;r')U=o = 0, 



(3) 



where r = {rii,z}. The solution of Eq. Q with the 
boundary condition Q is given by 

G(r,iV;r') = Go(r„ - A^; 0) [Go(z - z', A^; 0) 

- Gq{z + z',N;0)]. (4) 

It was argued in Ref. [23| that for an incompress- 
ible polymer melt in the presence of a neutral surface 
one should impose the reflecting (Neumann) boundary 
condition 



a,G(z,A^;z')U=o =0 



(5) 



on the Green's function of single polymer chains. The 
solution of Eq. Q with the boundary condition lO is 



n{z,N) = N, 



(11) 



respectively. The distance z in Eq. (|10|l is measured in 
units of Rg. The monomer density of one chain in the 
presence of a hard wall does not depend on the distance to 
the wall in the case of the reflecting boundary condition. 
The multiphcation of n{z, N) in Eqs. l(Tn|) and ((TT|l with 
the number of chains per volume n/V gives the monomer 
density of a mixture of independent chains. The necessity 
of change of the boundary condition in the polymer melt 
will be discussed at the end of Section ITlDI 



B. Collective description of the polymer mixture in 

bulk 

In the analytical approach to the description of dense 
polymer systems due to Leibler ^ the random phase ap- 
proximation is used to derive the Ginzburg-Landau type 
functional of the diblock copolymer melt as a functional 
of the order parameter. The collective description of 
concentrated polymer systems due to Edwards ,20j gives 
the physical quantities under interest as series in pow- 
ers of the effective potentials. These series are closely 
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related to those in the theory of polymer solutions [23 
with the main difference that the bare interaction poten- 
tials between the monomers are replaced by the effective 
ones (see below). The diagrammatic way of introduc- 
tion the collective description in Ref. |2J| enables one 
to go beyond the random phase approximation and es- 
tablishes the connection between Leibler's and Edwards' 
approaches. 

We now will consider the collective description of the 
polymer mixture consisting of A and B polymers in terms 
of effective potentials following Ref. |2J|, where this ap- 
proach was developed for copolymer melt. The elastic 
part of the Edwards free energy of ua polymers of type 
A and ns polymers of type B chains confined to a volume 
V is given by 



3 

2P 



UA+riB 

E 



N 



ds 



ds 



(12) 



where (s) parametrizes the configuration of mth poly- 
mer as a function of the position along the chain s. The 
interaction part of the free energy (in units of ksT) of 
the blend can be written using the microscopic monomer 
densities of both polymers 

riA+riB „N 

Pb{v) = E / ds5{v-v^{s)) (13) 



in the form 



where 



-Pint = i j d^ri j dV2Pa(ri)F„/3(ri -r2)p/5(r2), (14) 
ere 

V^,{v^^v,)^^y\^ ^ + ^)5(3)(ri-r2) (15) 

(a,/3 = A, B) is the interaction matrix of monomer- 
monomer interactions, and x is connected to the Flory- 
Huggins parameter. The sum convention over repeated 
indices is implied in Eq. I|14|) and henceforth. 

Let us now start with the computation of the average 
concentration of one of the polymers 



(Pa(r)) 



/ Dri{s)pa{Y) exp {-F^\ - F-, 



intj 



/ Dn{s) exp(-Fel-Fint) 



(16) 



using the collective description of the polymer blend in 
terms of effective potentials. The average monomer den- 
sity can be written after introducing a two-component 
field <i>Q(r) in the equivalent form as follows 

(Po(r)) 

SDn{s)jD<^>{v)Y{^,,5{<^p - p^)<i>„(r)e-^--^- 



(17) 



The insertion of the Fourier transformation of the infinite 
product of (5-functions 



\{6{<^p{v')-pp{v'))^ j DQ{v)e^' 

r',/3 



Q-(*-p) 



into Eq. 117|l . and replacement the order of integrations 

<a Q Q Q Q Q 
_i_ + _i i- + —i i + ... 



FIG. 1: Examples of graphs associated with the expression 

over the fields $(r) and Q{r) with the average over poly- 
mer configurations ri(s), yields the average over polymer 
configurations in the form 



J L>$(r)$„(r)e--^'" 



(18) 



where Q ■ p stands for J d^r Qa{r)pa{r), and the brack- 
ets (■..)g means the average over conformations of ideal 
polymer chains according to 



(19) 



To perform the average over polymer configurations we 
expand the first exponent in expression H19|) in Taylor 
series. The mean value (I19II decomposes as products 
of averages over single polymer chains, which have the 
structure 



dVi... / d3rfeQ„(ri)...Q„(rfe)(p„(ri)...p„(rfe)), (20) 



where k = 0,1, and a = A, B. According to Ref. |2J| 
it is convenient to associate expression (I20II with a graph 
containing k wavy legs, which are associated with Qa{j^i)- 
An example of graphs with k = 1, 2, 3 is shown in Fig. ^ 
The continuous lines are associated with the propagator 
(0) for a polymer blend in bulk and l@J for a polymer 
blend in the presence of a hard wall, respectively. Conse- 
quently, the series ((T^ is associated with a product of ua 
lines for A polymers and n b lines for B polymers contain- 
ing an arbitrary number of wavy legs in each line. Note 
that below we will consider the mean value H19|) in the 
thermodynamic limit ua 00, Ub 00, V 00 un- 
der the condition that the monomer densities computed 
using the one-polymer Green's function are constant: 

(p„(r))o = iV„n„/V = pa. (21) 

The corresponding density-density correlator reads 



(Pa(r2)Pa(ri))o = -^PaSaain - r2). 
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It is convenient to introduce the diagonal matrix Sap 
such that the Fourier transforms of the diagonal ele- 
ments, Saa{k), are the bulk structure factor of the ath 
component, which are given by 



(22) 



with g{y) — 2/?/^[exp(— y) + y—l] being the Debye func- 
tion. 

In order to carry out integrations over the two- 
component field Q(r) in treating a blend in bulk 24] one 
performs a partial summation of the series H19() (the lat- 
ter can be carried out only in the thermodynamic limit) 
by taking into account only the lines with one and two 
insertions (wavy legs in Fig. ^ in one polymer line. As 
a result one obtains the expression 



exp 



1 



• (a> -(/,)„)- -Q • 5 • Q 



(23) 



where (pq)q is the average monomer density The 
integrations over Q for a polymer blend in the bulk is 
easily performed in Fourier space and result in 



exp 



(24) 



where (5$(r) = $(r) — {p)q. The expression obtained 
after performing integrations over Q{r) can be written as 
exp(— with H{S^} being the Ginzburg-Landau 
functional — 

According to Ref. the functional integration over 
(5$ in Eq. ((T7|l yields the monomer density {pa{r)) as 
a series, which can be associated with Feynman graphs 
similar to those in the theory of polymer solutions in good 
solvent (see, for example, Ref. Q) with the difference 
that the bare potentials are replaced by the effective ones 

^ {V~^ + S)-\ (25) 

The lowest-order corrections to the monomer density are 
depicted in Fig.|21 The external lines in these graphs are 
associated with the expression 



(y + S-')-'S-\ (26) 

which can be written in the form V~^V'^^ . The con- 
tinuous lines denote the bare bulk one-polymer Green's 
functions 



C. Screened effective potentials in an 
incompressible polymer blend 

We now will consider in more details the properties of 
the effective potential (|25|l . We remind the reader that 
the quantities 5", V, and V^^ in Eq. are matrices. 
The elements of the matrix are explicitly given by 



vfAk) - n{k)[-v + 2VxSBB + x^SBB], 



vilik) 



VfBik) 



VfAik) - 

7^(fc)[-F- 



-n{k)[v ^ 

2VxSaa 



(27) 



FIG. 2: Examples of graphs contributing to the monomer 
concentration: graphs a and h are first order and c is second 
order in V^^ . After renormalization the continuous line is 
associated with the efltective propagator. Graph d with only 
one insertion of V^^ is identically zero after renormalization 
of internal lines. 



where for the sake of simplicity we have introduced the 
notation 



7^(fc) = [-I-VSaa-VSbb 

+ 2VxSaaSbb + X^SaaSbbV^- 



(28) 



The behavior of the effective potentials in polymer blends 
was studied in Refs. |23,l2al23- following we will 

consider an incompressible and athermic polymer blend, 
which in the formalism under consideration is described 
in the limit V oo and x ~^ 0. The effective potentials 
(|27|l simplify in this limit to 



vfAk) = 



1 



Saa + Sbb 



Using the explicit expression of the structure factor H22|) 
we obtain for large polymer chains 



VfAk) 



12 pa/I'a + Pb/I'b 



(29) 



As it follows from Eq. (|29|l the expansion in powers of ef- 
fective potentials is in fact an expansion in inverse powers 
of the density. 

We now will consider in more details the properties of 
the external potentials H26() associated with external lines 
in graphs a, b, and c in Fig. |21 which are explicitly given 
by 



T/CXt 

^AA 



(fc) = -nik)[l + VSBBh 

vxi^(k) = n{k)[SAA{v + x)], 
vE'l(k) = n{k)[SBB{v + x)], 

^SbW = -n{k)[l + VSAA]- 

In the case of athermic polymer blends the following iden- 
tities hold 

vT2{k) - vr^ik) = 1, y-i(fc) - v^-2ik) = i- (30) 
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For incompressible and athermic polymer blends V^^' (fc) 
simplify to 



s 



BB 



Saa + Sbb 



Sa/ 



Saa + Sbb 



For large polymer chains we finally get 
Pa/I-a + Pb/I-b 



Pa/Ia 



PaII\ + Ps/li 



and similar for V^)^{k) and Vg^(fc). Note that in this 
limit the external lines are independent of the wave num- 
ber k. Therefore, in the real space the external potential 
are local and are given by the Dirac's 5-function in this 
limit. 



a 



Q 



FIG. 3: Vertex with two insertions generated by the second 
term of Eq. CT . 



(Pa(r2)Pa(ri))o 

= T^Pa / ds2 / dsi {6[Z2 - z{s2)]S[zi - z{si)] 

^ JQ Jo 

X (5[r2ji - r2ji(s2)](5[ri,„ - r2,i,(si)])j, 



= / exp(i/C|,(r2,,i - ri.i,))S'QQ(z2, 21, iV),(31) 

where the Laplace transform of the diagonal element of 
the structure factor is given by 



D. Collective description of the polymer mixture 
in the presence of a hard wall 

We now will consider the collective description of a 
polymer blend in the presence of a hard wall. In contrast 
to the collective description of polymer blend in the bulk 
outlined in Sec. Ill Bl we have to use now instead of the 
free propagator given by Eq. (0) a propagator fulfilling 
an appropriate boundary condition. In a theory based 
on the statistical-mechanical description of single poly- 
mer chains, the boundary conditions should coincide with 
those of single polymers, i.e., be the Dirichlct boundary 
condition (^J. Since the behavior of a polymer chain 
in solution and in a polymer melt in the presence of 
a wall may be quite different, one can expect that the 
one-polymer Green's functions in solution and in melt 
may obey different boundary conditions. A consistent 
statistical-mechanical theory of polymer melt should be 
able, in principle, to derive the boundary condition for 
one-polymer Green's function appropriate for melt. We 
will show here that using a partial summation of graphs 
it is possible to reformulate the description of polymer 
blend in terms of the effective one-polymer Green's func- 
tion. We will bring forward the arguments that the latter 
should obey the reflecting boundary conditions. 

In order to introduce the collective description in the 
presence of the wall, we perform the same steps as in the 
bulk and arrive at the expression H18() . and expand fur- 
ther (exp(— • p))^ in Taylor series as given in Eq. ((^OJl. 
In contrast to the bulk the continuous lines are associ- 
ated with bare one-polymer Green's functions obeying 
the Dirichlet boundary condition Q. The field $a(r) as 
well as Qa{r) are defined in the whole space as in the 
bulk formalism. The mean density obtained by using the 
adsorbing boundary conditions is given by Eq. IIUI) mul- 
tiplied with the factor n^/V. The computation of the 
density-density correlation function {pa{ri)pa{'r'2))o (no 
summation over a) for a polymer blend in the presence 
of a hard wall gives 



Saaiz2,Zi,k„,p) = -^[Go(z2~' Zi, p+x)-Go{z2 + Zi, p+x)], 

(32) 

and where the notation x — kf^a? has been introduced. 
The nondiagonal elements of the matrix Sq,^ are zero. 
Note that mean values of density products are zero if 
one of Zi is negative, so that the expression (|32|l applies 
only for positive zi and Z2. 

For a polymer blend in the presence of a wall the trans- 
lationally invariant part of the structure factor (|32() is 
defined only in the half space, so that the integration 
over Q, which requires the inversion of Sap, is not so 
straightforward. In this case we separate the density- 
density correlator in two parts according to 

(Pa(r2)Pa(ri))o = (Pa(r2)Pa(ri))ob + (Pa(r2)Pa(ri))o 

- (Pa(r2)Pa(ri))ob 

= (Pa(r2)pa(ri))ob + (Pa(r2)Pa(ri))o, , (33) 

and perform a partial summation by taking into account 
in every line only the first term in Eq. (|33|) . In proceeding 
in this way we fix the reference state to be that of the 
bulk far from the wall. The prize to pay is that the 2nd 
term in Eq. 1)33(1 has to be taken into account as a vertex 
with two insertions, which is shown in Fig. 13 

The summation over lines with one and two insertions 
in one line results exactly in the expression given by 
Eq. H23|) with the average density given now by Eq. UlUfl . 
The terms in the series (|19|l with more than two fields 
Q{v) along one line [and two fields corresponding to the 
2nd term in Eq. (|33|l ] can be obtained from Eq. H23|l 
and consequently Eq. (|24|l as derivatives with respect to 
(^^(r). To compute the concentration profile according 
to Eq. (|17|l one should perform integration over the field 
<I>(r). While after integrations over Q the series (|19|l de- 
pends on (5$(r) = <I>(r) — (p), the interaction part of the 
free energy H14|l has the form _Fint = ^$ • F • $. Rewriting 
the latter in terms of (5$ yields 



1 



-S^ + d^-V ■ {p) . 



(34) 
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where the hnear term in 5$ has the same form as inter- 
action with an external field in the formalism of <i>^ the- 
ory. For an incompressible and athermic polymer blend, 
which we consider in the present work, the linear term 
vanishes. 

Similar to the consideration in bulk the expression ob- 
tained after performing integrations over Q{r) can be 
written as exp(— iJ{(5"I>}) with H{S^} being a Ginzburg- 
Landau Hamiltonian including the surface terms. How- 
ever, in contrast to the effective surface Hamiltonian used 
in many studies the corresponding terms are not lo- 
calized at the surface only . The integration over the 
field <I>(r) can be performed in the same way as for bulk. 

The collective description developed above is based on 
the concept of the effective potential, which takes into 
account the screening of monomer-monomer interactions 
in a melt. However, the effect of the wall is taken into 
account as in the case of diluted polymers via the Dirich- 
let boundary condition for one-polymer Green's function 
and leads to an inhomogeneous monomer density for dis- 
tances up to the gyration radius. However, in a melt 
the density is expected to be rather homogeneous at dis- 
tances z < Rg. This is the result of the interplay of the 
interaction with the wall and the incompressibility of the 
polymer melt. While in the polymer solution (which is a 
liquid and as such is incompressible) the entropic repul- 
sion with the wall favor the presence of solvent molecules 
at the wall. In the case of the melt this is not anymore 
the case, because the place of monomers being repulsed 
from the wall, will be occupied by monomers belonging to 
another polymer which at that moment are not or less re- 
pulsed from the wall. Due to this the melt density similar 
to the total density of the solution will not tend to zero in 
approaching the surface. We expect that the effect of the 
wall on the behavior in the polymer melt can be formu- 
lated in terms of the renormalized one-polymer Green's 
function, which should guarantee the uniformity of the 
density, and according to this should obey a boundary 
condition, which is different from the Dirichlet bound- 
ary condition. We now show that, indeed, the partial 
summation of graphs including insertions into continuous 
lines enables one to formulate the description of a poly- 
mer melt in terms of the effective one-polymer Green's 
function. We will consider for simplicity the renormal- 
ization of the bare one-polymer Green's functions in the 
expression for the concentration profile 

N °° °° 

= naj ds J dz'G{z',z,Na,- s) J dz"G{z, z" , s),{35) 

° 

where the bare one-polymer propagator obeys the Dirich- 
let boundary condition. The graphs b and c in Fig. [21 
contribute to the bare Green's functions in (|35|l . Using 
property H30|l and expressing the external potentials 
as 1 + we can divide this contribution into two parts. 
The first one is given by graphs b and c, in which the ex- 
ternal line associated with 1 and which renormalizes the 



one-polymer propagators in (|35|l . The second part, with 



z z z z 

E = -: r + — + 



FIG. 4; The lowest-order graphs contributing to the self- 
energy. The continuous lines are associated with the effective 
one-polymer propagator Gr.a- 

the external line associated with V^^^*, together with the 
graph a describes the fluctuations corrections to the con- 
centration profile. The renormalization to the first order 
can be extended to higher orders, with the result that 
the bare continuous lines will be replaced by the effective 
ones associated with the effective one-polymer Green's 
function. This procedure corresponds to reduction of the 
whole set of graphs to the skeleton graphs, i.e. the graphs 
without insertions into internal lines. The only exception 
are the graphs b and c in Fig. |2 which are due to the 
recasting of V^^. The renormalization of one-polymer 
graphs due to insertions into the internal lines can be 
represented using the Dyson equation 



G-^ = G-^ - E, 



(36) 



where E is the self-energy, which takes into account inser- 
tions along the chain. Note that G, etc., in Eq. 136|l are 
matrices with respect to spatial coordinates. Examples 
of graphs contributing to E are given in Fig. 2| 

As a result of the partial summation of graphs taking 
into account the insertions into internal lines according 
to Eq. H3t)|l the lowest-order contribution to the density 
profile (|35|l changes to 

oo 

{Pa{z)) = 71a y" ds J dz' Gr,a{z' , Z , Na - s) 



X J dz"Gr.a{z, z" , s). 




(37) 



The fluctuation corrections to Eq. H37|) are given by the 
skeleton graphs in Figs. |2 and [S] As a result of the par- 
tial summation the bare one-polymer propagators G are 
replaced by the effective ones Gr,a- 

Equation (|36|) with E given as an infinite set of graphs 
is the basis of the self-consistent computation of the effec- 
tive one-polymer Green's function in the polymer blend 
under presence of a hard wall. The solution of this equa- 
tion is a difficult task which goes beyond the scope of the 
present article. 

Fortunately, the form of Gr.a in polymer fiuid can be 
found from general arguments avoiding the direct solu- 
tion of Eq. H3()l) . According to the above discussion we 
expect that the density profile in an incompressible fluid 
in the presence of neutral wall will be uniform. On the 
another hand the density profile without taking into ac- 
count the fiuctuations is given by Eq. H37|) . As we have 
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A,z 

gI 



A,z 



VS'(r-r.) 



VAB(i--r2) 



FIG. 5: The Feynman diagrams giving the leading contribu- 
tion to the excess monomer concentration. 



shown in Sec. Ill Al the computation of the density using 
the one-polymer Green's function obeying the reflecting 
boundary condition gives a homogeneous density. Due 
to this we identify the effective one-polymer propagators 
Gr,a with that obeying the reflecting boundary condition. 
Deviations from Silberberg hypothesis in thin polymer 
fllms were studied recently in Refs. p^lsollsills^. 

FigureQlshows that the contributions to the self-energy 
does not reduce to an effective potential as it was as- 
sumed in Ref. in the approach based on the self- 
consistent field theory. While the first graph in Fig. 0| 
takes into account the monomer-monomer interactions 
along one polymer, the 2nd graph (and higher-order 
graphs) takes into account monomer-monomer interac- 
tions between different polymers. Since the graphs con- 
tributing to the self-energy take into account the many 
particle interactions characteristic for a melt, we ex- 
pect that the effective one-polymer Green's function Gr,a 
obeys the boundary conditions appropriate for an incom- 
pressible liquid, i.e., the refiecting boundary conditions. 
This conclusion is supported by the following argument. 
In polymer melts similar to semidilute polymer solutions 
the relevant quantity governing the properties of the sys- 
tem is the number of monomers between two subsequent 
cross-links along the chain, which for polymer melts is 
of order of unity, instead of the chain length N . Con- 
sequently, in the polymer melt the effect of the wall on 
the monomers, which are close to the wall, will be sim- 
ilar to that on solvent molecules in solution. However, 
the monomers are classical objects, which are described 
in the relaxational regime. For a single monomer (or 
a solvent molecule), which dynamics is described by the 
Langevin equation, the steady-state distribution function 
is given by Boltzmann distribution exp[— J7(z)/fcBT] and 
is therefore constant in the space between two walls. This 
makes clear that the Dirichlet boundary condition is ir- 
relevant in dense polymeric systems. 

It is well known that a polymer configuration corre- 
sponds to the trajectory of a quantum particle for imag- 
inary times. According to this the problem of bound- 
ary condition in polymer melts is expected to have its 
counterpart in quantum fluids in the presence of a neu- 
tral boundary. While the wave functions of single parti- 
cles obey the Dirichlet boundary condition at the wall, 
the density of the fluid is not required to be zero at the 
waff ml. 



III. COMPUTATION OF THE EXCESS 
MONOMER CONCENTRATION 

The skeleton graphs in Fig. [21 give the fluctuational 
part to the density proflle. The free end of the external 
line is associated with the argument of the monomer den- 
sity z (due to the symmetry along the wall the monomer 
density does not depend on r,,). The explicit calculation 
shows that the one-loop graphs, where the external line 
is located outside the loop (graphs b and c in Fig. [21) , are 
negligible for large N . The leading contribution is due to 
the graph a and the related graph, which describes the 
effect of B polymers on the concentration of A polymers. 
These graphs are shown in more details in Fig. [5l 

We now will consider the computation of the concen- 
tration of say the component A in the presence of a hard 
wall. We will assume that the statistical segment length 
of the polymer A is larger than that of the polymer B, 
I A > Ib, so that the polymer A is stiffer. The contri- 
bution to the excess concentration to the lowest order in 
powers of the effective potentials is given by graphs in 
Fig. [SI To conduct calculations it is convenient to con- 
sider the Laplace transform with respect to the contour 
length N. The analytical expression associated with the 
first graph in Fig. [SI is given by 



d^,, / dq VfAqf,+q^) 



. 1 



'.^2^-2zV^/a ^ 2gaVpT^sin((jz)e~^^'^/° 



P'^ {p + x){p + X + (p-a^Y 



(38) 



where p is Laplace conjugate to N and x = qf^d^ . The an- 
alytical expression of the 2nd graph in Fig. [3 is obtained 
from Eq. 138() using the replacements 

yTk - ^Tb. yfA - yfs. pa - pb. na - nb. 

Note that the factor —1 is due to the fact that V and 
appear with the sign minus in the exponential of 
the statistical weight of polymer configurations. The fc^ 
dependence of the effective potentials (|29|l leads to the di- 
vergence of the integrals over the wave vector in Eq. H38(l 
at the upper limit of integration. However, the effective 
potentials acquire for finite V their bare values for large 
k, so that the integral converges at the upper limit of the 
integration. Therefore, for finite V the effective poten- 
tials are screened only for lengths larger than the local 
length 



Ir 



Pb/I 



1-1/2 



which is obtained from the explicit expressions of the ef- 
fective potentials H27I) . The derivation shows that this 
length is the same for both polymers. We expect that 
for finite V the polymer blend can be considered as an 
incompressible only for lengths larger than Ic- In order to 
simplify the integration over the wave vector in Eq. H38|) 
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we use the athermic and incompressible limit of the effec- 
tive potentials H29|l . but restrict the integration to wave 
vectors smaller than the cutoff value A ~ 

The inspection of Eq. 1)381) shows that it (and the ex- 
pression associated with graphs with the external line 
being outside the loop) contains a z independent contri- 
bution to the excess concentration of the density. The 
straightforward computation yields the renormalization 
of the bulk monomer concentration as 



PA = PA 1 + (1 - 2) 



/2 ;2 



47r {pAll + PbIa^ I 



(39) 

The factor 2 in Eq. (|39|l accounts for graphs similar to the 
graphs b and c in Fig. [21 but with the external lines being 
on the right side of the interaction line. Note that the 
mass divergences of the graphs b and c are omitted, 
that implies the regularization of expression (|16|) with 
respect to the mass divergences at the beginning. Equa- 
tion H39|l shows that even in the bulk the packing effects 
change the bare density of the constituents: the concen- 
tration of the stiffer polymer becomes smaller. Without 
incorporating the possibility for a local nematic order- 
ing, which is not taken into account in the model of a 
Gaussian polymer chain, polymers with larger statisti- 
cal segment length are expected to have smaller density. 
Note that the renormalization of the bulk composition 
is local, and the comparison of Eq. H39|) with the corre- 
sponding expression for ps shows that the total density of 
the blend does not change. Although the renormalization 
of the bulk composition given by Eq. (|39|l is somewhat 
unexpected, its necessity can be explained qualitatively 
as follows. The density of an incompressible liquid at 
given T and V is determined by interactions between the 
molecules, and cannot be chosen arbitrarily as in gas- 
like systems. Thus, in application of the coarse-grained 
model under consideration to polymer blend Eq. H39|) 
describes the renormalization of bare concentrations to- 
wards their concentrations in the polymer melt, which 
are determined by monomer-monomer interactions. 

The z dependent part of Eq. H38() gives the excess 
monomer concentration as a function of the distance to 
the wall. The integration over the wave vector yields the 
simple expression 



X4) 



Sa^p^TTZ 



_e-(22/a)Vp+a"A"(a 



e-(22/")v^(a - z^) 
zp 



2z 



\/p + a?A^ 



(40) 



where Taix) — dt ^ exp{—t) is the incomplete 
Gamma function, and the notation 



Vo 



1 PaPb 



12 N All {paII\ + PbIII? 



is introduced. To obtain the excess density one should 
add to expression l|in|) , which is associated with the first 
graph in Fig. the corresponding expression associated 
with the 2nd graph in Fig. El 

We will first compute the excess concentration of the 
stiffer (A) polymer at the surface 5pa{z = 0). To that 
end we put z = in Eq. (|40|l . take into account the 
second graph in Fig. O and perform the inverse Laplace 
transform. For large N we obtain the result 



5pa[z = 0) 



3 paPbI\1% 



47r2 {paII+PbI\? 



^ ln(a|A2iVB) - \n{a\K^NA) 



.(41) 



If both polymers have the same gyration radius Rg 
aVN Eq. l|lT|l simplifies to 



PaPbI\1'\ 



1 



1 



(42) 

The excess concentration at the wall for polymer blend 
differing only in degrees of polymerization is derived from 
Eq. (gTl) as 



5pa{z = 0) 



3 A paPb 
47r2 P (pa + Pb) 



In- 



Nb 
Na 




FIG. 6: (Color online) Concentration profile of A polymers as 
a function of the distance to the surface for different values of 
Nb, and U = 1.5, Ib = 1, A"^ = 1.55, pa ^ Pb = 0.5. The 
continuous line: Na = Nb = 10''; dashes: Nb ~ 5x10^; dots: 
A^fl — 2x 10^. The inset shows the concentration profile in the 
vicinity of the surface as a function of the distance measured 
in units of Ia- 



The latter shows that the shorter polymers are present 
in excess at the wall. Notice that the excess concentration 
depends logarithmically on the number of segments N . 
The contribution to the excess concentration at z = 
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associated with graphs b and c in Fig. [21 reads 
1 PA 



PbI\? 



Due to the factor N^^^^ the latter vanishes for large N. 
Note that for conformationally asymmetric polymers of 
the same gyration radius the sign of Eq. H43|) is oppo- 
site to that of Eq. (|42|) . The increase of 6pa{z — 0) 
with N agrees qualitatively with the results of numerical 
simulations and calculations using the integral equation 
theory 4J. 

To compute dpA{z) for arbitrary z one should perform 
the inverse Laplace transform of Eq. 140|l . Since it cannot 
be performed analytically, we have used a numerical rou- 
tine (Durbin) for inverse Laplace transform in Mathemat- 
ica. The results of the numerical calculation of the excess 
concentration of stiffer polymers Spa{z) for different val- 
ues of the degrees of polymerization of more flexible poly- 
mer are shown in Fig. It shows that the increase of Nb 
results in an increase of the excess concentration of the A 
polymer. For Nb < Na the concentration of A polymers 
is still in excess in the vicinity of the wall, but becomes 
lower than in the bulk for intermediate distances, i.e., 
the B polymers are in excess at these distances. These 
results are in agreement with numerical simulations and 
computations using the integral equation theory . Fig- 
ure [7| shows the result of the computation of the excess 
concentration of the shorter polymers in a polymer blend 
consisting of chemically identical polymers, which differ 
only in their degrees of polymerization. Figure shows 
that shorter polymers are present in excess in the vicinity 
of the wall. This finding is in qualitative agreement with 
the result predicted in Ref. [s^ and observed in Refs. Q 
and IsEImIsI. 
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FIG. 7: (Color online) Concentration profile of A polymers as 
a function of the distance to the surface for different values of 
Nb, and U = Ib = L5, Na = 10*, A"' = 1.55, pa = Pb ^ 
0.5. The continuous line: Nb =2x10"; dashes: A^s = 5x10*. 



sideration is compatible with the excess of the solvent at 
the wall in a polymer solution. The latter corresponds 
to the limit, when the polymerization degree of shorter 
polymers tends to unity. However, to describe this limit 
one has to take into account the higher-order terms in 
the perturbation series for the concentration profile. 

Note that the both cases we have considered above 
{I A ^ Ib, RgA = RgB and Ia ^ Ib, Na ^ Nb) follow 
from the general formula (|40|l . 

We now will give a qualitative explanation of the dif- 
ferent behavior of polymers in the blend under the influ- 
ence of a hard wall. A single polymer in a dilute solu- 
tion obeys the Dirichlet boundary condition. As a conse- 
quence of the boundary condition the number of config- 
urations available to the polymer chain lowers with the 
decrease of the distance to the wall. The latter results 
in an entropic repulsion of the polymer from the wall, 
and is responsible for the vanishing of the density at the 
wall. According to this the solvent molecules are favored 
in the vicinity of the wall with respect to the polymer 
monomers. A simple calculation using the distribution 
function obeying the Dirichlet boundary condition shows 
that the force acting on the free end of the polymer at 
a given distance to the wall is controlled by the gyration 
radius of the polymer Rg — a\fN . 

A completely different behavior takes place in the case 
of incompressible polymer melts, where the entropic re- 
pulsion from the wall is balanced by the melt pressure 
with the consequence that the density is uniform. How- 
ever, there is a difference in the behavior of the polymers 
in the vicinity of the wall for melt composed of different 
polymers. We consider first a polymer blend composed 
of polymers which differ only in degrees of polymeriza- 
tion. In a layer with the thickness equal to the gyration 
radius of larger polymers, the larger polymer experiences 
the entropic force from the wall while the shorter poly- 
mer does not. Due to this the larger polymer increases 
its distance to the wall, which will be occupied by shorter 
ones, in order that the total density will remain constant. 
The asymmetry in the behavior of polymers in the vicin- 
ity of the wall appears even in a polymer melt composed 
of identical polymers. According to the above argument 
the monomers of a polymer coil, which has contacts with 
the wall, are disfavored with respect to the ends of poly- 
mer coils which do not have contacts with the wall. Due 
to this the polymer ends are expected to be present in 
excess in the vicinity of the wall. The effect of the distri- 
bution of polymer ends on the surface tension was studied 
in Ref. '3»| . A quantitative study of the distribution of 
polymer ends using the self-consistent field theory was 
performed in Ref. |39l |. 

For polymers with different statistical segment lengths, 
but the same gyration radius the difference in the behav- 
ior in the vicinity of the wall can be explained qualita- 
tively as follows. The monomer density of a polymer coil 



is given by pc = N / ~ a 
sity of a coil is ps = pcRg = 



"/Rg, while the surface den- 
■ a~^. Therefore, the surface 



The excess of shorter polymers in the case under con- density ps of the stiffer polymer is smaller. It is likely to 
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expect that the repulsive effect of the wall on the coil is 
proportional to Ps ■ According to this the repulsive effect 
of the wall is stronger for more flexible polymers. This 
is the reason that the monomers of stiffer polymers will 
be favored in the vicinity of the wall. The surface enrich- 
ment 6pA is expected to be proportional to the differences 
of surface densities, i.e., SpA ~ P^~pt' which agrees with 
our quantitative result (|^ . According to this qualitative 
consideration the difference in surface densities pf — p^ 
is a drive for the conformational asymmetry. Since the 
monomers within the layer of thickness Rg are affected 
by the wall, we expect that the excess concentration will 
depend on Rg. However, the logarithmic dependence on 
Rg in Eq. (|42|) is difficult to derive using only the hand 
wavy arguments. 

Note that in the above computation of the excess con- 
centration SpA^z) we have taken into account the lowest- 
order correction in the series in powers of effective po- 
tentials. The effective potentials according to Eq. H29|) 
are inversely proportional to the density, so that the per- 
turbation expansion in powers of effective potentials is a 
series in inverse powers of the density. However, since the 
polymer melt has a fixed density, the inverse density is 
not a small parameter. The magnitude of the first-order 
correction can be controlled by considering polymers hav- 
ing the same gyration radius and small differences in Ia 
and Ib, or polymers with small differences in Na and 
Nb for I A = Ib- However, it is not clear without explicit 
computations, if the 2nd order term is smaller than the 
1st order one under the above conditions. From the gen- 
eral point of view one would expect the following bounds 
on the total effect of the perturbation series. As already 
mentioned above for polymers differing only in degrees 
of polymerization the effect of the whole perturbation 
series should recover in the limit Na <Si Nb the behavior 
in polymer solutions, where the polymer concentration 
will tend to zero in approaching the surface. For poly- 
mers differing in flexibility the concentration of the stiffer 
polymer at the wall cannot exceed the total density of the 



polymer blend in bulk. In other words the concentration 
of the more flexible polymer cannot be negative. This 
determines the upper limit of applicability of our results 
given by Eqs. (UIJ) and (O- 



IV. CONCLUSIONS 

To summarize, we have generalized the Edwards' col- 
lective description of dense polymer systems in terms of 
effective potentials to polymer blends in the presence of 
a surface. Using this formalism we have studied an in- 
compressible athermic polymer blend of conformationally 
asymmetric polymers, which differ in statistical segment 
lengths, in the presence of a hard wall. We have com- 
puted the excess concentrations of constituents to the 
first order in powers of effective potentials. We have 
found that stiffer polymers are in excess in the vicinity 
of the surface, and that the concentration excess at the 
surface depends logarithmically on the degrees of poly- 
merization. For polymer blends differing only in degrees 
of polymerization the shorter polymers are in excess at 
the wall. Our results are in agreement with numerical re- 
sults available in the literature. The present method can 
be applied in a straightforward way to study the behav- 
ior of polymer blends and copolymer melt in the presence 
of selective surfaces, to study the dimensions of polymer 
molecules in the melt, the distribution of polymer ends, 
etc. 
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